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lewis Research Center 
Cleveland, Ohio 


SUMMARY 

Open-circuit voltages measured in silicon solar 
cells made from 0.01 D-ca material are 150 mV lower 
than voltages calculated from simple diffusion theory 
and cannot be explained by poor diffusion lengths or 
surface leakage currents. An analytical study was 
made to determine whether high doping effects , which 
increase the intrinsic carrier concentration, could 
account for the low observed voltages and to determine 
the limits on voltage and efficiency imposed by high 
doping effects. The results indicate that the observ- 
ed variation of voltage with base resistivity is pre- 
dicted by these effects. A marimiaa efficiency of 19^6 
(AMI) and a voltage of .7 volts were calculated for 
0.1 3 -an cells assessing an optimum diffused layer im- 
purity profile. 


INTRODUCTION 

There has been considerable interest in raising 
the efficiency of silicon solar cells to near its lim- 
it. Analyses ( 1 , 2 ) have been made to define what 
gains can be expected and also what areas of research 
are meet likely to yield the necessary improvements 
(or to identify the principal obstacles that must be 
overcome ) . 

It has been suggested that the maximum practical 
silicon solar cell air mass zero efficiency of 1 8 % 
would occur in 0.1 ohm- cm material. Realization of 
this efficiency depends heavily on obtaining an open- 
circuit voltage close to 0.7 volt calculated from 
simple diffusion theory. The best present-day per- 
formance of low resistivity cells, however, falls sub- 
stantially short if these predictions. Msasured effi- 
ciency and voltage are typically about 10^ and 0.6l 
volts, respectively. It is believed that improvement 
of the voltage is the key factor leading to high effi- 
ciency. However, the difference between the achieved 
and predicted voltage performance is one area for 
which a satisfactory explanation is not available. 

It has been widely hypothesized that the low ob- 
served performance of the 0.1 ohm- cm cell is a result 
of undesirable carrier recombination processes occur- 
ring in the base and depletion regions. However, meas- 
urements show that the short-circuit current and the 
base diffusion length are not primarily responsible 
for the low cell performance. A study of low resis- 
tivity cells by Soclof and lies (3) not only support 
this , but also strongly indicate that other mechanisms 
which influence junction injection efficiency must be 
considered. One such mechanism is the narrowing of 
the silicon bandgap by high dopant concentrations. 

The theoretical work of Kane ( 4 ) and Morgan (5) 
demonstrates that donor or acceptor concentrations in 
excess of 10^° atoms/cm^ alter the conduction band and 
broaden 4 ne impurity band densities of states. This 
work has been applied to silicon material by several 
investigators ( 6 - 10 ). Their results suggest that high 
impurity concentrations can indeed change the bandgap. 
Additional theoretical and experimental studies are 
cited in (10). 

An important consequence of bandgap narrowing for 
p-n junction device operation is the increase of the 


intrinsic carrier concentration JU . Van Overstraeten 
et al. (11) have incorporated in detail the impurity 
level-dependent in the basic semiconductor device 
transport equations. Similar equations modi' - ' »d to 
include bandgap narrowing have been used to t alyze 
and reconcile the long-standing discrepancy between 
the predicted and measured current gain in bipolar 
transistors (6, 12-16). 

The alteration of the bandgap and the intrinsic 
concentration by high iapurity concentration effects, 
which will be referred to as the ffeavy Doping Effect 
( HDE) in this paper, may also have significant effects 
on silicon solar cells. 


The purpose of this paper is to make a preliminary 
analytical assessment of the significance of the HDE on 
the open-circuit voltage of silicon solar cells. Volt- 
age calculations based on diffusion theory (such as 
used in ( 17 )) are modified to include the first order 
features of the HDE. Comparisons of the open-circuit 
voltage measured for cells of various base resistivi- 
ties are made with those calculated using the diffusion 
model with and without the HDE. 

THEORETICAL BACKGROUND 
Open- Circuit Voltage 

The open-circuit voltage (V ) is calculated using 


V 
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where the diode saturation current (I 0 ) with the HDE, 
given by 
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is the sun of the diffused and base region components 
( Ip ■ I^j + I ofe )- The " d " and " b " subscripts refer 
to the diffused and base regions , respectively. AH 
calculations are made for an n + -p device and assume 
ohmic contacts, uniform n+ region Impurity level, a 
0.25 m junction depth, a 3 00 urn cell thickness, and 
an operating temperature of 300 K. For the purpose of 
this paper, the effect of bandgap narrowing on the 
short-circuit current (I ) is neglected. Values of 
I aC used in the calculations are typical of the various 
cell resistivities. The mobilities (^, m^) and impur- 
ity concentrations (N_., N— ) corresponding to material 
resistivities were selected from the literature ( 18 ). 

A diffused region diffusion length (L^) of about 1.0 
is assumed. However, the base region values are 
based on typical experimental data. The geometry fac- 
tors and are given in detail in ( 17 ). However, 
for the case considered here they reduce to coth W/L 
(i. e. infinite surface recombination velocity) contacts 
are ass\saed. 
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For those calculations which include the HDE, the 
values of corresponding to a net Impurity concen- 
tration are selected from the calculated data of (11). 
These data are shown in figure 1. It is seen that for 
net impurity levels below to* 8 / cm3, n has the classi- 
cal value (n^ c ) of I.lxl0 10 /cm3. More Importantly . 
hcvever , n^ increases rapidly for levels above 10*9. 
Typical solar cells have diffused region impurity con- 
centrations in excess of 10*9/cnr, thus the variation 
of n^ must be taken into account. The implication of 
figure 1, suggested in (4,5) is that the conduction 
and valence band edges are no longer well-defined at 
high impurity levels. Hence the bandg&p is narrowed 
and n^ Increased, compared to a lightly doped, slight- 
ly extrinsic crystal. 


Calculations which do not include the HDE are re- 
ferred to as the simple diffusion theory. For this 
case equations (l) and (2) are also used; however, the 
intrinsic carrier concentration is independent of 
impurity level and equal to the classical value (n tc ) 
of l.lxlO^/cm 3 . 


Effective Inpurlty Profile 

An effective impurity concentration is defined by 


Hft = 


"ia^ 




(3) 


1)1?) i a the actual donor concentration at seme location, 
x in the material and the ratio, n id ( N id)/ n ic* i* a 
heavy doping factor based on the data or figure 1. 

Using equation, (3) the diffused region saturation cur- 
rent component (Iod) C9jri written in the form given 
by 


n| c (kT)(Md) 
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Note that for 


"id < 10 ' 
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Because 1^ becomes large when is small, the dif- 

fused region properties can dominate the voltages cal- 
culated from equation (l). 

Four layer Model 

Calculations of the voltage, current, and effi- 
ciency were made for several p-base resistivities us- 
ing the four layer model (17,19). These calculations 
include those features that are currently believed 
necessary to optimize cell performance, e.g. , X_ = 

0.1 m 9 3% reflection losses, a p + back contact, and 
the optimum n + impurity profile (an actual gauss ian 
profile with a 1x10*9/ cnP surface concentration). Dif- 
fused layer and p-base lifetimes were based on experi- 
mental data (3). A diffused layer lifetime of 3x10"^ 
sec corresponds to an impurity level of lxlO^/caP. 

The p-base diffusion lengths are 303, 206, 201, 100, 
and 39 ^an for 10, 1, 0.5, 0.1, and 0.01 ft-cm resis- 
tivities, respectively. The rear surface recombina- 
tion velocity (SRV) is 10? cm/sec. Data are calculat- 
ed for two values of front SRV, 1 q 3 and 10? cm^/sec. 
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RESULTS ADD DISCUSSION 
HDE on Impurity Profiles 

It is clear from equation (3) that as the actual 
doping level Increases, the effective doping concentra- 
tion, as modified by the HDE, will decrease. Thus sub- 
stantial alterations to observed diffusion profiles 
would be expected. Figure 2 indicates the influence of 
the HDE on a gauss ian profile having several different 
aurfa"^ concentrations. For concentrations above about 
lCr°/cit3 > the effective impurity profile is retrograde. 
The slope of the retrograde portion increases as the 
surface carrier concentration increases. Thus it would 
appear that in these cases, an opposing drift field is 
set up whi^h prevents flow of minority carriers to the 
junction. This explanation is an alternative to the 
Head layer* model of the cell. The retrograde profile 
can be eliminated by reducing the surface concentration 
to about 10*9 7 cm3, a gauss ian profile results with an 
effective surface concentration of 3x10^/ an3, Gaussian 
profiles are obtained for surface concentration less 
than or equal to lO*9/cm3. The electric field (11) as- 
sociated with this gauss ian effective profile is such 
that the minority carriers throughout the region are 
directed toward the collecting junction, thus aiding 
carrier collections. It appears that optimum cell per- 
formance will result for the maximum doping concentra- 
tion at which the effective profile h:.s no retrograde 
sections. 

HDE on Open- Circuit Voltage 

Results of two types of calculations are to be 
presented. First, those of voltage using equations (l) 
and (2) assuming a uniformly doped n + diffused layer. 
Second, an overall cell performance calculation using 
the four layer model with the gauss ian profile and other 
parameters described above. 

Figure 3 summarizes the highest measured voltages 
for 10, 1.0, 0.1 and 0*01 ohm-an silicon cells. Also 
shown as the solid line are the "simple theory" voltages 
calculated with no HDE and where 1^ is insignificant 
canpared to I 0 v for all base resistivities. As seen, a 
discrepancy exists between the calculated and measured 
voltages which becomes more severe with decreasing base 
resistivity. The discrepancies are about 75 mV and 
200 mV for 0.1 and 0.01 ohm- cm cells, respectively. 

Uniform Doping Calculations . Using the values for 
in figure 1 and equations (1) and (2), the voltages 
are recalculated for the 0.1 and 0.01 ohm- cm cases for 
various values of uniform diffused region impurity lev- 
el. Since the p-base properties for each resistivity 
were not varied, the changes in the total I Q and the 
voltage were due only to the different diffused layer 
impurity levels. These calculated data are ihown in 
figure 4. For comparison, the simple thc^ 1 results 
are also shown for a 0.1 ohm-cm cell. It is seen that 
the simple-theory V Qc rises, then saturates with in- 
creasing impurity level. This result is expected since 
it reflects the diminishing influence of 1^ 

Calculations which include the HDE contrast sharply 
with the simple theory. The voltage is seen to rise 
initially but, rather than saturate, peaks at about 
0.66 volts, for a concentration of 1x10*9 atoms/cm . 

The V then decreases reflecting the influence of o 
large values i.e., the HDE. Note, that at 2xl0 £ ‘°/cm^, 
which corresponds to the present day diffused layer im- 
purity level, the calculated voltage of nearly .60 volts 
is in agreement with the experimental data of figure 3‘ 

A trend, similar to the 0.1 ohm-cm case, occurs 
for the 0.01 ohm-cm cell. Although the voltage again 
peaks at an inpurity level of about 1x10*9 atans/cnp^ 
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its magnitude is only 0.60 volts. Interestingly, a 
voltage of 0. i6 volts is predicted for an impurity 
level of 2x10 . This voltage alao compares well with 

the measured value shown in figure 3. 

A aumnary of the experimental and peak calculated 
voltages is shown in figure 5. The open circles and 
the solid line are the data of figure 3* The dotted 
line represents the peak voltages calculated with the 
HDE for the uniform diffused layer impurity concentra- 
tions. The significance here is that in addition to 
predicting the experimentally observed trend of v oc 
with base resistivity, the HDE theory also suggests 
that voltage improvements appear possible by reducing 
the doping level in the diffused region. 

Four layer Model Calculations . A more detailed 
analysis of cell performance was made using the four 
layer model. Calculations of efficiency, voltage and 
current as a function of base resistivity are shown 
in figure 6. It is seen that for the 1CP SRV case a 
peak efficiency of nearly 19$ AMO and a V QC of 0.7 V 
occurs for 0.1 G-cm material. These values are in 
agreement with previous estimates (l). They also sug- 
gest that HDE need not reduce cell performance provided 
that the optimum diffused layer doping profile can be 
achieved . 

For the 10^ SRV case, these calculations yield a 
lower peak efficiency of 16.2$ and a lower voltage of 
about 0.6 V. The base resistivity for peak efficiency 
shifts to 1 D-cm material. This decrease in performance 
is most probably related to the value of lifetime cho- 
sen for the diffused layer. Use of bulk material life- 
time values would have improved the performance in this 
case. There is a strong sensitivity of these calcula- 
tions to the shape of the profile, the value of SRV and 
the lifetime in the layer. Of these parameters, least 
is known about the lifetime in the diffused layer. This 
parameter is critical to separation of the various 
models for predicting solar cell performance. Thus a 
research effort should be undertaken to determine this 
lifetime and thus provide necessary input data for 
theoretical calculations. 

SUMMARY AND CONCLUSIONS 

This paper describes a preliminary analytical 
assessment of the influence of heavy doping on the 
open- circuit voltage performance of silicon solar 
cells. The modification of simple diffusion theory, 
guided by published calculations of the intrinsic car- 
rier concentration as a function of impurity concentra- 
tion, can explain the experimental voltage trends in 
low resistivity cells. Furthermore, the calculations 
suggest that voltage improvements can be expected, 
provided that the diffused layer impurity concentration 
and profile are optimized. 

An assessment of the maximum practical efficiency 
for 300 ixn thick cells in the 0.01 to 10 ft-cm range 
was also made using the four-layer model. These calcu- 
lations assume a p + back, 0.1 junction depth, 3$ 
reflection losses . front surface recombination veloci- 
ties of 103 or 105, a nd the optimized diffused layer 
profile. A maximum AMD efficiency of about 19$ occurs 
for 0.1 D-cm material. This supports earlier predic- 
tions and suggests that heavy doping effects should 
not be deleterious to cell efficiency providing the 


diffused layer impurity profile is optimized. Fc * the 
lifetime values used, cell efficiency was sensitive to 
surface recombination velocity, decreasing to 16$ for 
SRV * 10- . Cell performance is critically dependent on 
the profile, the value of SRV and the lifetime. Be- 
cause lifetime is the least well known, a research ef- 
fort to obtain the necessary data is required. 
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N d -N a (ATOMS CM' 3 ) 

Figure 1. - Effect of dopant concentration on 
intrinsic concentration. 
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DISTANCE INTO DIFFUSED LAYER 
(ARBITRARY UNITS) 

Figure 2. - Effective net impurity profile in diffused 
layer of 0. 1 Q-cm solar cell (for actual Gaussian 
profiles of various surface concentration). 
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Figure 3. - Experimental and predicted open-circuit 
voltage dependence on base doping levels. 
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Figure 4. - influence of diffused layer heavy doping 
effects on open -circuit voltage. 
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Figure 5. - Experimental and predicted open-circuit voltage 
dependence on base doping level. 
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BASE RESISTIVITY (fl-CM) 

Figure 6. - Calculated cell performance as a function of base 
resistivity (using the four layer model and two values of 
front surface recombination velocity). 
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